arXiv: 1506.02714v2 [cond-mat.soft] 20Jun2015 


Free volume under shear 


Moumita Maiti/ H. A. Vinutha,^’^ Srikanth Sastry,^’^ and Claus Heussinger^ 

^Institute for Theoretical Physics, Georg-August University of Gottingen, 
Friedrich-Hund Platz 1, 37077 Gottingen, Germany 
^ Jawaharlal Nehru Genter for Advanced Scientific Research, 

Jakkur Gampus, Bengaluru - 560064, India 
^TIFR Genter for Interdisciplinary Sciences, 

21 Brundavan Golony, Narsingi, 500075 Hyderabad, India 
(Dated: June 23, 2015) 

Abstract 

Using an athermal quasistatic simulation protocol, we study the distribution of free volumes 
in sheared hard-particle packings close to, but below, the random-close packing threshold. We 
show that under shear, and independent of volume fraction, the free volumes develop features 
similar to close-packed systems - particles self-organize in a manner as to mimick the isotropically 
jammed state. We compare athermally sheared packings with thermalized packings and show 
that thermalization leads to an erasure of these structural features. The temporal evolution, in 
particular the opening-up and the closing of free-volume patches is associated with the single¬ 
particle dynamics, showing a crossover from ballistic to diffusive behavior. 
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I. INTRODUCTION 


The transition from a flowing liquid to a jammed solid state has been a subject of research 
in a wide range of systems, including granular matter, colloidal suspensions, and diverse 
types of glass formers, and in the context of gelation. A particular context, in which the 
geometry associated with such a transition is important, has been that of the jamming to 
unjamming transition in athermal particle packings, often modeled using sphere packings, 
with or without the presence of external driving. The relationship between the jamming 
phenomenology and thermal systems, either the rheology of driven thermal ^sterns, or to 
the glass transition in undriven systems, are intensely investigated currently [ll-l3|. 

Disordered assemblies of spheres undergo a jamming transition at a packing fraction of 
~ 0.64, at which the pressure for hard sphere packings diverges. At this random close 
packing (RCP) density, individual hard spheres are entirely constrained by their neighbours 
and have no space to move around. The network of such spheres in contact (“backbone”) 
spans the system and may thus support external load. 

The geometrical properties of sphere packings and their influence on the mechanical re¬ 


sponse is a complex problem with many different 


acets. One key observable is the connec- 
. Wyart [8| discussed the consequences of 


tivity of the network which is isostatic at RCP jd-G] 
opening of contacts, as well as the distribution of gaps, i.e. distances between non-contacting 
particles. Atkinson et ah js] have studied the structure of rattlers, the particles that are 
not constrained by enough contacts. Schroder-Turk et ah [10| have observed a signature of 
RCP in the shape of Voronoi cells of the particles. 

Packings at RCP also have special mechanical properties. The bulk modulus, for example, 
is hnite for soft-sphere packings at RCP, while it generically vanishes for spring networks of 
equivalent connectivity [ll, 12| Apparently, the particles in a packing can organize in order 
to resist compressive forces, in a way that is not possible for spring networks. 

In the vicinity of random close packing, spheres have finite free space to move, and 
free space vanishes approaching RCP. For thermal systems this implies slow dynamics 
approching RCP. At the same time, dynamics in sheared systems at zero temperature 
speeds up by approaching RCP. Indeed, particles in shear flow are found to move ever faster. 


the less space they have 


13| . It is an interesting but essentially unsolved question how such 


non-trivial dynamics arises and couples to other properties of the system, in particular to 
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FIG. 1: Probability distribution of free volumes (2D): (a) (p = 0.825 for decompressed (unsheared) 
configurations and for different elapsed strains; (b) Different volume fractions in the steady state 
and compared with the distribution for nearly jammed configurations, (c) Data for the 3D system. 


the singular behavior of the correlation length or the rheological coefficients 
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181. Here, 


we are interested in the analysis of the free space available for particle movement, and how 
that may be affected by externally imposed deformation, and in turn how such changes in 
geometry influence particle flow. 

The divergence of pressure in a hard particle system is a consequence of the vanishing 


Q. 


of free space. Pressure is the ratio of free surface to free volume [711, the latter vanishes at 
jamming. In the fluid phase, the free volume distribution is peaked around the mean free 
volume, and at jamming the distribution can be expected to have a delta function peak at 
zero free volume for the backbone spheres. Indeed close to jamming, for isotropic, thermal- 


ized, ensembles of configurations at any given density, the free volume distribution 


19| has a 


strong peak around the mean free volume. A surprising observation, however, is that there is 
a power law tail in the free volume distribution in addition, which appears to be a signature 
of being close to the jamming point. In the present work, we address the question of how 
the free volume distributions of sphere packings close to the jamming point get modihed 
under shear deformation, both in the athermal and thermal cases, and what we may deduce 
regarding the organization of spheres in response to external shear. 


II. SIMULATION DETAILS 

We study two- (2D) and three-dimensional (3D) systems of N soft frictionless spheres. 
To avoid crystallization, we use two different particle sizes; ^ spheres have diameter d and 
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Y have diameter lAd. The particle density is quantihed via the fraction of volume (area in 
2D, although we use ‘volume’ in the 2D case as well, when no confusion is caused) 0 that is 
occupied by the particles. 

Particles interact via elastic repulsive forces. Two particles repel each other with a 
harmonic potential energy. 


^el = -6(1 - 

CLi-i 


Tij dij 


( 1 ) 


where is the distance between the two particles. The cut-off = {di + dj)/2 is set by 
the diameters of the two interacting spheres. 

To implement the shear a quasistatic simulation protocol is used j^. B |. In the qua¬ 
sistatic limit (small-strain rates, 7 ^ 0 ) the average overlap <5 = (1 — ^), vanishes and the 
particles effectively behave as hard-spheres (equivalent to e —?■ 00 ). In recent work 18, Q it 
was shown that the statistics of particle velocities obtained from the quasi-static simulations 
is identical with the small-strain rate (Newtonian fluid) limit of fully dynamic molecular dy¬ 
namics simulations. Thus, our quasistatic results should be considered representative for 
the Newtonian flow-regime of dense suspensions close to jamming. 

The jamming density under shear is at 0c = 0.647 (0c = 0.842 in 2D). We probe densities 
below this limit in the range 0 = 0.61 — 0.64 (0 = 0.825.. .0.840 in 2d). In this range of 
densities cooperative effects set in and the correlation length increases roughly by a factor 
of ten 13]. 


III. RESULTS 

Motivated by recent results on the free volume distribution in isotropically jammed pack¬ 
ings |l^, we ask how steady shear affects these distributions. The free volume of a particle is 
the volume (area in 2D) that is available for the center of the particle with all other particles 
are hxed in space. 

Starting with a jammed conhguration at 0c, we decompress the conhguration to the 
target density at 0 < 0c by reducing the diameters of all particles, di —)■ adi, where a 
is the scaling factor for the particle diameters. In this way, all particles acquire a certain 
amount of free volume Vf (or free surface Sf in the two dimensional case) with a probability 
distribution, P{vf), that is highly peaked at Vf ~ a'^ (see Fig. [1] (a) for the 2d system, blue 
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line/downward triangles). In this figure we also highlight the distribution of free volume for 
the nea rly j ammed configuration (before decompression), which is a power-law (black open 
circles) |l9 |. 

The “decompressed” configuration serves as starting point to our shear simulation. Fig. [T] 
(a) indicates the evolution of the probability distribution as a function of the elapsed strain. 
After only small amounts of strain, all particles retain their finite free volume, but as strain¬ 
ing goes on the number of those particles reduce and a peak at zero free volume starts 
developing (not shown in the figure). At the same time, a broad tail develops. In the 
steady-state, thus, the form of the distribution is rather different than in the decompressed 
configuration. Instead of a narrow, peaked distribution we obtain a strong peak at zero free 
volume with a broad power law tail. 

Fig. [1] (b) illustrates that, close to jamming, the power-law tail is independent of density 
0, which is also identical in the exponent to the isotropically jammed state. The same 
phenomenon is observed also for the three-dimensional system, as shown in Fig. [T] (c), 
illustrating that the outcome is not dependent on dimensionality. Apparently, during shear 
particles self-organize in such a way as to mimick the isotropically jammed state - eyen 
though the system is not jammed and flows as a Newtonian liquid. 


It is worthwhile comparing the situa 


iments performed at low densities 
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ion of steady shear with that of cyclic shear. Exper- 


25| have shown that after sufficient shear cycles the 


system self-organizes into an “absorbing” state, where particles do not interact anymore. 
Their trajectories during a shear cycle are strictly reversible. This indicates that particles 
generate free volume around themselves in such a way that during their shear-induced os¬ 
cillations no interactions with other particles occur. Under steady-shear the organization 
apparently is opposite, and free volume is not generated but destroyed (Fig. [2]). The aver¬ 
age free volume under shear is roughly two orders of magnitude smaller compared to the 
decompressed configuration at the same density. 

Another illustrative comparison is with thermal systems. Recall that the system we study 
is inherently athermal. Particles only move because they are driven by shear. To generate 
a thermalized system at the same density we again use the decompressed state as starting 
configuration, but now run short Metropolis Monte-Carlo simulations. The resulting free 
volume distribution is presented in Fig. |3]and compared to the sheared state. 

It is clear that thermalization acts rather differently than shearing and does not produce 
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FIG. 2: Average free volume as a function of volume fraction. Comparison of configurations under 
shear and decompressed states. Two different system sizes are compared to exclude finite-size 
effects. 



FIG. 3: Comparison of free-volume distributions from jammed, sheared and thermalized configu¬ 
rations (top: 2D, bottom; 3D). 


a power-law tail in the free volume distribution (in agreement with Ref. 
interesting also from the point of view of shear-induced effective temperatures 
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This is 


29|, which 


rely on the assumption that driving by shear is in some sense equivalent to thermal driving 
but at an “effective” temperature. At least for the observable under consideration - the 
free-volume distribution - this equivalence does not seem to hold here. 

Another difference between the sheared and thermal systems is the distribution of the 
total free volume among the particles. In the thermal system (and also under cyclic shear) 
the free volume is distributed basically among all particles. In the sheared system, only a 
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certain fraction of particles have a non-zero free volume. In the density range we studied, we 
hnd a fraction of 6% to 11% of such hnite free volume particles. The majority of particles 
belong to a backbone of locally jammed particles with zero free volume. The full free-volume 
distribution therefore consists of a power-law tail for particles with hnite free volume plus a 
delta-function peak for the backbone particles. As the density is increased, more and more 
hnite free volume particles are incorporated into the backbone, such that at (pc the backbone 
becomes globally jammed and the system solidihes. 

A structural self-organization via shear how is also visible in the pair-correlation function. 
The power law singularity as well as the splitting of the second peak are commonly taken 
as the signature of the jammed state. Fig. 0] shows that sheared conhgurations exhibit 
a similar power-law singularity with the same exponent as in the jammed state, whereas 
thermalization destroys the divergence. Also the splitting of the second peak in steady- 
shear is as sharp as in the jammed state but it is smeared out under thermalization. We 
have also computed the radial distribution function only for the hnite free volume particles 
and hnd it to be similar as for particles with zero free volume. This is consistent with the 
hndings of Atkinson et al. 9|] who report a certain amount of spatial correlations in the 
rattlers. 

Solidihcation is known to occur when the backbone is exactly isostatic. We can extract 
the approach to isostaticity also from the free volumes. To this end we shrink the particles 
(in the steady-state shear conhgurations) by a very small factor /. This generates a hnite 
free volume v = 0{f'^) for the backbone particles. The average number of facets of these 
free volume patches then dehnes the connectivity of the backbone. 

Fig. [5] shows that the backbone network is hypostatic (as expected in the huid regime 
below jamming) and thus not mechanically stable. Isostaticity is approached as ziso — z ^ 
{pc — PY, with an exponent x close to one. Being hypostatic, there are a number of zero 
frequency modes along which particles can move without cost in energy. This motion is, 
for example, visible as a short-time ballistic regime in the mean-square displacement j^. 
Similarly, this motion is rehected in the temporal evolution of the free volume of the hnite 
free volume particles. The locations of such particles represent holes in the network of 
backbone particles and how these holes close or open-up rehects how particles move in the 
network. 

We calculate the strain needed to close a hole. This represents a measure for the “life- 
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FIG. 4: Pair correlation function g{r) for jammed packing, as well as sheared and thermalized 
configurations (color code as in Fig. [3|). Sheared configurations display the same characteristic 
power-law slope as the jammed configuration. Thermalization destroys this power-law. Inset shows 
that sheared configurations have a split-second peak as seen in jammed configurations, whereas for 
thermalized configurations, the sharp features are smeared out. 



FIG. 5: Effective connectivity in the backbone of zero-free volume particles - measured via the 
number of facets of the free volume after tiny amount of decompression. 

time” of finite free volume for a particle that resides in this hole. The probability distribution 
of strains to close the holes is displayed in Fig. [6l The striking feature is its power-law 
tail, which is a consequence of the broad distribution of hole sizes. Holes can close after 
arbitrarily small strains. On large strain-scales, the tail is cut-off at strains 7 c that decrease 
with increasing the volume-fraction towards (pc- In fact, the cut-off 7 c is comparable in 
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FIG. 6: The distribution of strain values needed for closing a hole. The arrows signify the strain 
7msd; at which the particle dynamics (as measured by the mean-square displacement, msd) crosses 
over from ballistic to diffusive (taken from Ref, [l^ i. 


magnitude to the strain 7msd, at which the mean-squared displacement has a crossover from 
ballistic to diffusive behavior (extracted from Fig. 1 in Ref. 13|). The strains for three 
different (j) are indicated in Fig. [6] by arrows. This connects the evolution of free-volume 
patches - holes in the network of backbone particles - with the single-particle dynamics. 
The cross-over to diffusion signals the onset of the decorrelation of the floppy, zero-energy 
modes, and this decorrelation is associated by the opening and closing of these holes. 


IV. CONCLUSION 

In conclusion, we have shown that unjammed athermal hard-particle systems under shear 
self-organize such that structural properties resemble those of jammed packings. In par¬ 
ticular, we have reported the probability distribution of free-volumes as well as the pair- 
correlation function. 

The free-volume distribution in (nearly-)jammed packings has previously been shown to 
display a power-law tail [l^ . We show here that this tail is also present in sheared systems 
in a hnite interval of densities below the jamming threshold. By way of contrast, thermalized 
conhgurations in this density range generically show a peaked and rather narrow distribution 
of free volumes. 

Apparently, shearing and thermalization act rather differently what regards these lo- 
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cal structural properties. Shearing drives the system “towards” the jamming point, while 
thermalization drives the system away from it. 

The average free volume in the sheared conhgurations is roughly two orders of magnitude 
smaller than in the thermalized conhgurations. Moreover free volume is heterogeneously 
distributed among only a few hnite free volume particles, with the majority of particles 
belonging to a locally jammed backbone with zero free volume. 

Our results raise the question about the relevance of hnite free volume particles for the 
how properties close to jamming. Indeed, a power-law in the free-volume distribution means 
that there are particles with arbitrarily-small amounts of extra space. As a consequence, 
they can easily be integrated into the backbone after only inhnitesimal amounts of strain. 
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